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Introduction 



The tetrahedron equation is the three dimensional generalization of the Yang-Baxter 
equation which guarantees the existence of commuting transfer matrices. The 
importance of Yang-Baxter equations for modern mathematics and for mathematical 
physics is well known. However, the nature of the tetrahedron equation is much less 
understood, this mainly because it is a much more complicated equation. 

Given the physical interest to understand the nature of the singularities which 
give rise to 3-D phase transitions, any effort which gets us closer to analytic results 
for 3-D statistical systems seems worthwhile. What has been achieved recently, is to 
construct large classes of 3-D solvable models with Z^-spin variables and to streamline 
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the otherwise complicated formalism. The goal of an analytic calculation of partition 
functions and order parameters is not yet close by. The only available result from Baxter 
[2E] does not lend itself to generalizations. 

The first solution of the tetrahedron equation was obtained in 1980 by A. 
Zamolodchikov Pfl 0] and then generalized by R. Baxter and V. Bazhanov [Hj (ZBB- 
model) and others [3]. These models have Z^-spin variables and solve the IRC 
(Interaction Round a Cube) version of the tetrahedron equation. Later also the solution 
of the dual equation, the vertex tetrahedron equation, was obtained [Sj, generalizing 
several vertex solutions known previously Ej- Here we shall consider only vertex- 
type solutions, which are usually denoted by R, the symmetry will include aZji. In 
general these R-matrices obey the so called "simple modified tetrahedron equation" 
which recently has been investigated in [Oj. The modified tetrahedron equation (MTE) 
allows us to obtain the ordinary tetrahedron equation for composite weights or vertices. 
In the IRC formulation this has been shown in [TU1 ITT] , while the most simple vertex 
case was considered in |12j . 

In this paper we shall introduce a new convenient theta-function parameterization 
of general R operators. This parameterization will allow us to define fused weights 91, 
which are partition functions of open cubes of size M 3 , and which obey a MTE. In 
special cases the 91 solve an ordinary tetrahedron equation. 

The vertex matrix 91 e End (C 37Vm2 ) is parameterized in the terms of iV-th roots 
of theta-functions on the Jacobian of a genus g = (M — l) 2 compact algebraic curve 
T g . The divisors of three meromorphic functions on T g play the role of the spectral 
parameters for 91. An additional parameter of 91 is an arbitrary v G Jac(T 9 ). The 
tetrahedron equation for 9^ holds due to M 4 simple modified tetrahedron equations. In 
the case when M = 1 and therefore T g = S2, the solution of the simple tetrahedron 
equation of |5] is reproduced. 

This paper is organized as follows. In section 1 we recall the vertex formulation 
of the 3-D integrable ZBB model and sketch the derivation of the matrix operator 
Hijk from a current conservation principle and Z-invariance. It satisfies the MTE 
and can be parameterized by quadrangle line-sections. In section 2 we introduce the 
parameterization of Ry/. in terms of theta functions. The Fermat relation and the 
Hirota equations are written as Fay identities. In section 3 we show that a theta 
function parameterization allows a compact formulation of the fusion of many R^ 
to the Boltzmann weight 91 of a whole open cube which satisfies a MTE. In section 4 
we first consider the special case of vanishing Jacobi transforms, in which 91 satisfies 
a simple TE, then we discuss the rational case and the relation to Chessboard models. 
Finally, section 5 summarizes our conclusions. 

1. The R-matrix and its parameterization 

In this first section we give a short summary of basic previous results which also serves 
for fixing the notation. 
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1.1. The H-matrix of the vertex ZBB model 

We start recalling the vertex formulation of the ZBB- model jH]. We consider a 3- 
dimensional lattice with the elementary cell defined by three non-coplanar vectors 
e 1; e 2 , e 3 and general vertices 

n = rii ei + n 2 e 2 + n 3 e 3 , m, n 2) n 3 G Z. (1) 

We label the directed link along starting from n by (j, n). On these links there are 
spin variables crj n which take values in Z^r. The partition function is defined by 

Z = J^J (°"l,n, 0"2,n+e 2 ; °"3,n| R l^l^+eu 0"2,n! C"3,n+e 3 ) • (2) 

W * 

where R is an operator (which in the ZBB model is independent of n) mapping the 
initial three spin variables to the three final ones, so that 

<^25 = K^IRK,^) (3) 

is a N 3 x N 3 matrix independent of n. 

For the vertex ZBB-model, Q can be expressed as a kind of cross ratio of four 
cyclic functions W p (n). Introduce a two component vector p = (x,y) which is restricted 
to the Fermat curve 

x N + y N = 1. (4) 

Then define the function W p (n) by 

n 

W p (0) = 1, W p {n) = ] [ _ V - for n > 0. (5) 

1 (J x 

where 

q = e 2 ^' N (6) 

is the primitive N—th root of unity Because of the Fermat curve restriction, W p (n) is 
cyclic in n : 

W p {n + N) = W p {n). 

Now R = H(pi,p 2 ,P3,P4) is defined by the following matrix function depending on four 
Fermat points Pi,P2,P3,P4 

°°*+^ 2 +«3 1 W Pz {<j' 2 - at) W P4 (a 2 - a[) ' 1 ] 

where x-coordinates of four Fermat curve points in (JJJ) are identically related by 

So the matrix elements R^\% 2 2 %1 depend on three complex numbers. These correspond 
to Zamolodchikov's spherical angles in the IRC- formulation of the ZBB-model The 
structure of the indices of the matrix (JJJ) allows one to consider R as the operator acting 
in the tensor product of three vector spaces 

V = , R G End (V ® V ® V) . (9) 
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It is conventional to enumerate naturally the components of the tensor product of several 
vector spaces, so that (JTj) are the matrix elements of R = R123. Of course, Ri 23 acts 
trivially on the all other vector spaces if one considers V (g,A for some arbitrary A. 

(|7ji is known as the R-matrix of the Zamolodchikov-Bazhanov-Baxter model, see 
[Sj. The proof that (|ZJ) satisfies the Tetrahedron Equation 

E ^S(p (2) ) ^Sr(p (3) ) ^S> (4) ) 

ji-je 

ii-ie 

is rather tedious 0. In (fTU|) the arguments (j = 1, ... ,4) stand for four Fermat 
curve points {P\\pi\pi\pA > ) each. These 16 points depend on five independent 
parameters expressible in terms of spherical angles, see Note that here on the left 
and right hand side the same pV> appear. This will not be the case in the generalizations 
which will be discussed soon. 

Baxter and Forrester (THj have studied whether this model describes phase 
transitions. They used variational and numerical methods and found strong evidence 
that for the parameter values for which (fTUj) is satisfied, the ZBB-model is just at 
criticality ^H]- So, in order to get a chance to describe phase transitions while staying 
integrable (recall that also for the 2D Potts model this is a problem), one should enlarge 
the framework and define more general Boltzmann weights and introduce less restrictive 
Tetrahedron equations. Less restrictive and still powerful generalized equations can be 
used, as has been shown by Mangazeev and Stroganov [TUj: they introduced Modified 
Tetrahedron Equations which guarantee commuting layer-to-layer transfer matrices. 
Further work along this line has been done in fH\ \1'2\. 

1.2. H-matrix satisfying the Modified Tetrahedron Equation 

Since in the above-mentioned work [3] the proof that particular Boltzmann weights 
satisfy a particular MTE has been rather tedious, here we shall follow the approach 
introduced in [21] in which there is no need for an explicit check of the MTE. The 
Boltzmann weights are constructed from "physical" principles which guarantee the 
validity of the MTE and nevertheless leave much freedom to obtain a broad class of 
integrable 3D- models. We give a short summary of the argument. 

One starts with an oriented 3-D basic lattice. The dynamic variables living on the 
links i of this lattice are taken to be elements iti,Wi G K>i an ultralocal Weyl algebra 
W = to.j at the primitive N-th root of unity: UiWj = q Sij WjUi , (q as in eq.fjHJ)), 
which generalize the Z n spin variables of (j2J) and (JHJ) . The Weyl elements are represented 
by standard N x N raising respectively diagonal matrices. The N-th powers of the Weyl 
variables are centers of the algebra and so are scalar variables. 

The main object constructed is an invertible canonical mapping IZyk in the space 
of a triple Weyl algebra. Tl^k operates at the vertices of the 3D lattice, mapping the 
three Weyl elements on the "incoming" links onto those on the "outgoing" links, see 
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Figure n The construction of TZ^k is based on two postulates, a Kirchhoff-like current 



i 




Figure 1. The six links of the basic lattice intersecting in the vertex A, intersected by 
auxiliary planes (shaded) in two different positions: first passing through roi, tt>2, tt>3 
and second through tn^, tv' 2 , rOg. The second position is obtained from the first by 
moving the auxiliary plane parallel through the vertex A. The Weyl variables, elements 
of tt>i, tnj live on the links of the basic lattice. 1Z maps the left auxiliary triangle onto 
the upper right one. 



Y 1Z 12 3 Y 




Figure 2. The canonical invertible mapping IZ123 shown in the auxiliary planes 
passing through the incoming (left) and outgoing (right) dynamic variables which are 
elements of ttii resp. ro^. The directed lines X, Y, Z are the intersections of the three 
planes forming the vertex A of Figure ^ Their sections are labeled by the line-section 
parameters 61, , . . . , d^. Note that the choice of the orientation of the lines is not 
unique. The orientation chosen here corresponds to the numbering l|44|) and l|45|l of 
the fused vertex considered in the Section 

conservation and a Baxter Z-invariance, and gives a unique explicit result: a canonical 
and invertible rational mapping operator. Since q is a root of unity, TZijk decomposes 
into a matrix conjugation Ry&, and a purely functional mapping 1Z\A which acts on the 
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scalar parameters (the Weyl centers). So, for any rational function $ on tv: 

7Z 123 o $ = R 123 (7^g o $)R- 1 3 . (11) 

It turns out that the matrix TUjk has the form (J2J) where the four Fermat curve 
parameters, again constrained by JHJ), are rational functions of the scalar Weyl center 
parameters. 

Next consider an auxiliary plane which cuts the three incoming links near a vertex, 
and a second auxiliary plane cutting through the outgoing links, see Figure H We take 
the six Weyl dynamic variables to sit on the six intersection points of the auxiliary 
planes. 7^123 can be regarded as the mapping of the ingoing auxiliary plane to the 
parallel shifted outgoing auxiliary plane. 

Now consider the vertices of the basic lattice to be formed as the intersection points 
of three sets of non-parallel planes. The three planes which form the vertex A of Figure^ 
intersect the auxiliary planes in the lines X, Y, Z shown in Figure El In Figure d these 
intersection lines are the sides of the shaded triangles. Seen from the moving auxiliary 
plane, 72-123 shifts the line X through the vertex with index 1 or Y through the vertex 
2 etc. We attach variables b±, b 2 , . . . , d 3 to each section of the lines X, Y, Z as shown 
in Figure 121 

It is convenient to parameterize the two scalar variables associated with the 
incoming dynamic Weyl variable tni (corresponding to u± and w± in usual notation) 
by the ratios c 2 /c 3 and d 3 /d 2 . Analogously, e.g. those for ro 2 are defined as /b' 2 N 
and d' 2 N /di etc. Details of the rule to parameterize the scalar variables in terms of 
"line-section" variables bx,...,ds etc. are explained in jH]. However, these will not 
be essential here, since one of the aims of this paper is to introduce and use another 
parameterization. Just observe in Figure El that 71^23 changes only three of the line- 
section parameters: b 2 , c 2 , d 2 . From the explicit form of the canonical operator 7^123 
(see |9 ) one finds that the functional mapping is rational in the N-th powers of 
the line sections: 

b»<$d? + b^C^d^ + K?b£c»d» 



b 



fN 



r ' N 



r N d N 
c 2 u 2 

K N h N r NjN , K N h N r NjN , „jV K N uN „N jN 

1 3 1 2 ' 3 2 3 1 ' 1 3 ^3 *"2 1 

Wbjd? ; 



uN NrfN , uN NjN , NuN N jN 

j fN 2 c l "3 ' °l c l a 2 i ft 3 °1 c 2 a l / 1 n \ 

"2 - TW-Jf • W 

u 2 c 2 

Here K\, k 2 , k 3 are fixed parameters ("coupling constants") of the mapping 7^i 2 3. One 
can show |Hj that in the line-section parameterization the three independent Fermat 
curve parameters which determine R123 according to (|7 ji .(|5 ji are 

b 2 c 3 K 2 bid 2 61C3 

xx = —r — ; x 2 = — ; x 3 = ■ 13 

K163C2 6 2 c i \/Q b 2 c 2 

If we define a partition function in analogy to (J2J), the matrix elements of Rijk take 
the role of generalized (because generically they will be complex) Boltzmann weights 
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of integrable 3D lattice models of statistical mechanics. Despite their non-positivity we 
shall just call these matrix elements "Boltzmann weights". 

Via the physical assumptions made in constructing IZijk , the validity of the TE is 
already built in. Simply considering two different sequences of Z-invariance shifts in a 
geometric figure formed by four intersecting straight lines ("quadrangle"), one concludes 
that (see e.g. P) 

7^123 ' 7^-145 ' 7^-246 ' 7^-356 ~ ^356 • 7^-246 ' 7^145 • 7^123 , (14) 

i.e. that TZijk satisfies the Tetrahedron equation. Inserting (JTTj) into (fT4*|) and choosing 
various phases of N-th roots (the Fermat points (fT3*j) involve b±, . . ., whereas the (fT2|) 
relate , . . .), leads to the MTE for the matrix operator R.; 



R 



123 



7^-123 R-145 



7^1237^145 R-24f> 



,(/) 
'145 



'M23 'M45 /c 246 u 



R356 j 



R 



356 



7^356 ° R 246 



7^3567^246 ° R 145 



'^356 '^246 'M45 ° ^123 



(15) 



Via the Fermat points each R^ depends on several scalar variables, see e.g. (|13J) . In 
(j!5|) the scalar variables which appear in the matrices TUjk are to be transformed by the 



( f) 

functional transformations 1Z)A as indicated. Let us write shorthand 



R« 

R( 4 ) 

R (6) ,v ,v 
lv — /v 356 /v 246 

Then (}20|) becomes 

R(D R(2) R (3) R (4) 



Rl23; 

npU) npif) TpU) n -R . 
7< -123 7< -145 7< -246 ^-356) 

'R'J^rR^AR Rl45! 



R( 2 ) 
R( 5 ) 
R( 7 ) 



7^-123 R-145; 



-!?(/) 7? (/)??(/) n"R ■ 

/■Cc>c;k /Co,,,; /-C^c; O rtl23, 



R( 3 ) 



7^1237^-145 R-246; 



-356 /v 246 /v 145 
7^-356 R 246; 



R 



356- 



(16) 



(17) 

C N . p is a scalar 



p# R( ? ) R( fi ) R^, 

where each R( J ) acts non-trivially in only three of the six spaces V 
density factor which comes in when passing from mappings to matrix equations. 

The parameters which determine the R^ are the corresponding Fermat curve 
coordinates. Taking into account the functional transformations in in terms of 
the line-section parameters one finds (for full details see [Oj): R^ = R(x^ ,x 2 \x$ ) 
with 



x 



(2) 



sn 



T (i) 

x 2 

T (2) 



r (7) 



;r 



(i) \ 

Q 



;r 



(2) 



X 



V 



2; 



(8) _(8) 



X 



3 

(8) 



-1/2 



/ 



/ &2C3 


K 2 b 1 c' 2 


&1C3 


K1&3C2 

a 2 c 2 


b' 2 d ^ 


aic 2 


Kia 3 ci 


a 2 c 


gV2 a5ci 


«2&3 


K A a\b 2 




n 2 a 3 b 2 


a}6{ 


gV2 a tt6 2 


a x b 2 




a Q b 2 


\ K 3 a 2 bi 


Mi 


gV2 ajfei 



The once or multiply transformed parameters like c'/, follow from the iteration 
of equations like (|12|) . Altogether, since there are eight matrices R^) appearing 
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in the MTEs, and as seen in Figure |2l each transformation changes three line- 
section parameters, we have 24 equations for 32 different line-section parameters (these 
parameters can seen in Table ^ below). This is a set of classical integrable equations 
which conveniently are written in Hirota form: 

b 2 ' N c N 2 d N 2 = b"<£d» + b N 2 c N ,d N z + K?b»cM; 

K 2 ^2 C 2 d 2 = &i b 3 c-l ' d 2 + n^b 2 c^d 1 + k x K^b^ c 2 d^; 
b^c^d 2 ' N = b?c»dg + fof cf d» + «f6f c^<; 

af cf < = of cf < + afcfV* + «f a^cf rf 2 /iV ; 

^4 ^2 ^1 1 — ^1 ^3 ^0 1 ~~ ^~ ^5 ^2 ^2 ^0 ~~ ^1 ^5 ^3 ^1 ^0 ' 

„N„NjllN _ „N^Nj iN i „N^NjN i „N„N„NjN. 

u> 2 C-y (1^ — (1 2 Cg li 2 i" C-q li^ "T Kg tig , 

af X d 2 /JV = a»b 2 ' N d 2 ' N + d%a?b 2 > N + ^d^af; 
n N ,a^ N d 2 ' N = ^b%afd 2 ' N + 4a^b 2 ' N d'[ N + ^a^&fdf; 
a^d 2 "' N = l$d»a? + a£b%d 2 ' N + K»a»b»df; 



n Vi N h IN IN _ h N uiN IN , h N N UN , IV iV JV, /JV 
1*2 "2 C 2 — "3 ^2 ~r "3 ^3 "2 ~r K3 CI3 63 02 



N ju ,iAT ■. inA7 in AT-, in AT t AT AT in AT, t AT f/JKf 



ifb 2 ' N cf = cfaf&T + a'? N b'? N c 2 N + K%a'? N b 2 ' N c" N - 



b'fcfdf = b N Q cf N df + b[ N cf N df N + K?bfc\ N df N ; 

K>2 ^1 ^1 ~~ ^1 ^2 ^1 H - ^3 ^1 ^2 ^1 ^3 ^0 ^2 ^1 7 

bfcfdf = c"bfdf N + b»c»df + K^d^cf. (19) 

The first three of these equations are just (fT2"j). defining TZ 123 , i.e. b' 2 N , c' 2 N , ef^ in 
terms of the unprimed &i, . . . , cf 3 . The first six equations together (e.g. express in the 
fourth eq. on the right hand side c' 2 N and d' 2 N from the first and third equations) define 
7^123 7^145' e ^ c - The complete expressions for (JTHJl and (fT9"|) can be found in 

Straightforward combination of the first twelve equations of (11 9 j) on one hand, and 
of the last twelve of equations of (fT§j) on the other hand (best done e.g. by Maple), 
shows that the functional mappings given in (|T9*|) automatically satisfy the functional 
TE: 

/v 123 /v 145 /v 246 /v 356 — /v 356 /v 246 /v 145 /v 123 l zu / 

where for the superposition of two operators acting on a function $ we use the notation 
((A ■ B) ■ <&) (^4 ■ (i3 • $)). Of course, the validity of (j^Uj) is a consequence of the 
physical rules used when constructing IZyk- In the line-section parameterization the 
relation between the first, second etc. lines in both (jl8j) and in (119(1 is not transparent. 
Introducing a new parameterization in the next subsection will make these relations 
simple and explicit. 
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2. Parameterization using concepts of algebraic geometry 



2.1. Theta functions 

It is well-known [TU EE H3 EE] that Hirota-type equations can be identically satisfied 
by a parameterization in terms of theta functions on an algebraic curve. We shall now 
introduce such a parameterization in order to write (|18j) and (jl9j) in a more systematic 
way. This will be useful also later to formulate fusion in a transparent manner. For the 
notations of algebraic geometry see e.g. [T3] . 

Let T g be an abstract generic algebraic curve of the genus g with u> being the 
canonical ^-dimensional vector of the homomorphic differentials. For any two points 
X, Y G T g let l£ : T] h-> Jac(r 9 ) be 



I 



Y 



(21) 



Y 



Let further E(X,Y) = — E(Y,X) be the prime form on T 2 , and 0(v) be the theta- 
function on Jac(r ff ). 

It is well known, the theta-functions on the Jacobian of an algebraic curve obey the 
Fay identity 



G(v) 6(v + 1 A B + Ig) = 9(v + l£) 9(v + 1% 



+ e(v + i£) e(v + i c D 



E(A, B) E(D,C) 
E(A,C) E(D,B) 
E(A, D) E(C, B) 



(22) 



E(A,C)E(D,B) ' 

which involves four points A,B,C,D G T g and a v e Jac(r 9 ). We shall show that in 
the parameterization to be introduced below, the Fermat relations become just Fay- 
identities. The Fay identity involves only cross ratios of prime forms, and these ratios 
have a simple expression in terms of non-singular odd characteristic theta functions: 



X X' 
Y Y' 



def E(X,Y)E{X r ,Y>) 



e 



e odd 



E(X, Y>) E(X', Y) Q eodd ( I*' ) 6 £odd ( II, ) ■ (23) 

For solving the 24 trilinear equations (|19p we shall need an identity with more arguments 
Q, X, Y, Y', Z, Z' G T g obtained by combining two Fay identities: 

e(v + iS)e(v + i? + i?)e(v + ]2 + i£) 



- e(v + ti + if) e(v + 1^) e(v + i w z + i Y Y 



- e(v + 15^ + 1?) e(v + + if) e(v + ti 



Z'\ 



+ e(v + \ Q X + i Y Y + if) e(v + 1?) e(v + ti 



Z'\ 



X Y 

Z Z' 

X Z 

Y Y' 

X Y 

Z Z' 



X Z' 
Y Y' 



. (24) 



Furthermore, since we will need the iV-th roots of theta-functions and prime forms, we 
also define e(X, Y) and 0(v) by 

e(X,Yf = & eodd (I*) ~ E(X,Y); 9(v) N = 0(v) . (25) 
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Since in the following we shall have to write many equations involving theta-functions, 
it is convenient to introduce special abbreviations. For Q, A, B x , B[,...eT g we define: 

n 

(A, B 1 + B 2 + ... + B n ) =e(v + l2 + £l§'); (A, B) = E(A, B) ; 

3=1 
n 

[A, Bl + B 2 + ... + B n ] =0( v + lj + ^lg'); 

3=1 

{A, B} = - g -V> e (A,flO/e(A5); {AS} = -e(A', S) / e(A, B) . (26) 

Note that the brackets ( , ) and [ , ] introduced here do not show explicitly the 
dependence on the variables v, Q, B[, ... B' n since these always come in the same form. 

We also introduce, using these notations: 

F (v; X, Y', Y, Z>, Z) t f (X)(Y, Z)(Z, Y) (Y, Z) (X, Z') (X, Y') 

- (X, Z) (Y) (Z, Y) (X, Z) (Y, Z') (X, Y') - (X, Y) (Y, Z) (Z) (X, Y) {¥', Z) (X, Z') 

+ (X,Y + Z) (Y) (Z) (X, Z) (X, Y) (¥', Z') , (27) 

so that the Double-Fay-identity pij) is 

F(v;X,Y',Y,Z',Z) = 0. (28) 

The dependence on Q is trivial since it appears only in the combination v + 1^ . So 
Q is not an independent variable. 

2.2. Re-parameterization of R 

Let us introduce the new parameterization of the matrix (J7J). As we illustrated in 
Figure in the auxiliary plane the mapping 1Z\2Z can be considered as a relative shift 
of three directed lines X, Y, Z with respect to each other. Now, for the given algebraic 
curve T g and v G C 9 , we introduce three pairs of points on T g : 

X', X, Y', Y, Z', Z G T g . (29) 
Another point Q G T g will just serve as a trivial normalization. Then let 

R = R( Pl ,p 2 ,p 3 ,p,) ^ R = R(v;X',X;r,Y;Z',Z) (30) 

with, using the shorthand notations (J2T?|) and pj = (xj,yj): 

_1 {X,Z} [X,Y][Y,Z] e(Z,Z')e(X,Y') [Z, Y] 0(v + g + if ) 

Xl q{Y',Z}[X,Y + Z][YY Vl e(X, Z) e(Y>, Z>) [X, Y + Z] [Y] ' 

_ {X',Z} [X] [Y,X + Z] e(Z,Z')e(X',Y) [Z, X] fl(v + g + if) 

X2 {Y, Z} [X, Z] [Y, X] ' V2 q e(X>, Z) e(Y, Z>) [X, Z] [Y, X] 

_1 {X,Z}[X] [Y,Z] e(Z, Z') e(X, Y) [Z] 9(v + ly + if) 



q {Y, Z} [X, Z\ [Y] ' tfi * e(X, Z) e(Y, Z>) [X,Z][Y] 

_1 {X',Z} [X,Y] [Y,X + Z] e(Z,Z')e(X',Y') [Z, X+Y] 0(v + ly + if 



g{r',Z}[A,r + Z][F,A]' e(X',Z)e(y',Z') [X, Y+Z] [y,X] 

(31) 
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Actually, see (j5J), for defining R 12 3 we don't need the jji themselves but only the ratios 



(33) 



y 3 = {Y,Z'} [Z] [X,Y + Z] y 4 _ {X,Y'} [Z,X + Y] [Y] y 3 _ {X,Z} [Z][Y,X] 
Vl q {X,Y} [X,Z][Z,Y] ' yi {X,Z} [Z,Y][Y,X] ' y 2 {X, Y} [Y] [Z, X] 

(32) 

from which e(Z, Z') and 6{y + Iy + if) drop out. 

Note that for this parameterization we used a generic algebraic curve and generic 
points on this curve, and a generic point on its Jacobian, all in order to parameterize 
just three independent complex numbers Xi,x 2 ,x 3 . In (|3~T]) all Xk, Vk are the periodical 
functions of v G Jac(lT 9 ) . 

The parameterization (|3*T|) is suggested by a few assumptions: First: The prime 

forms shall appear in the Xi only in the form of N-th roots of (|2*3*j) : 

r i l/W 

{X,Z} X Y 

{Y,Z} ~ Z Z' 

Second, considering (|T3^1. we demand that the line-section parameters b±, b 2 , b 3 , b' 2 
(sections of the line X in Figure |2J should be proportional to N-th roots of theta 
functions of the form [X, . . .] defined in (J2Ej). Analogously, the sections c\, . . . , c' 2 of 
the line Y are assumed to be proportional to [Y, . . .] . Finally, we consider that we want 
to use Fay identities to provide the Fermat relations and the Hirota equations. 

The merit of this parameterization will be seen in several places: when we consider 
the transformed mappings R( 2 \ R (6) , when we re- write eqs. ()19j) and when we 
construct composite weights in Sectional 

We now must verify that (|31|) . and its generalization to the other Fermat points 
in the MTE, give a consistent parameterization of the relevant equations (fTHj) and 
()19|) . We first check that (J31)) satisfies the Fermat relations 

x? + yf = I- (34) 

Indeed, these are true due to the Fay identity, which for A, B,C,D e T g we write as 

- (A, C) (D,B) 0(v) 9(v + li + Ig) + (A,B)(D,C) 6(v + 1^) 6(v + ig) 

+ (A,D)(C,B) 6(v + 1^) 8(v + Ig) = 0. (35) 

For j = 1 put in (JB5J) (A, B, C, D) -> (Y 1 , X, Z', Z ) and v -> v' = v + ig , giving 

(Z, X) (Y 1 , Z') (Y) (X, Y + Z) — (Z, Z>) (Y 1 , X) (Z, Y) 9(v' + if) 

- (Z',X){Y',Z)(X,Y)(Y,Z) = 0, 

for j = 2 put in (JSHJ) (A, B, C, D) -> (X', Y, Z', Z) and v -> v" = v + I J , giving 

(Z', Y) (X', Z) (Y, X) (X, Z)-(Z, Y) (X', Z') (X) (Y, X + Z) 

+ (z,z')(x',Y)(z,x)e(v" + i^) = o, 

for j = 3 put in © (A -B, C, D) —> (Z, X, Z', F) and v -> v+ = v + l Q z : 
(X, Z') (Y, Z) (X) (Y, Z) - (X, Z) (Y, Z') ( Y) (X, Z) 

+ (Z,Z')(X,Y)(Z)0(v + I« + l|') = 0, 
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or (X, Z')(Y, Z) 0(v+ + if) 0(v+ + if) - (X, Z)(Y, Z') 6(v+ + if) 6(v+ + if) 

+ (Z,Z')(X,F)e(v + )e(v + + lf + If) = 0. 
For j = 4 put in © (A B, C, D) ^ (Y', Z, X', Z'), v -> v* = v + ij + if . 

£.3. Line-section parameters and Hirota equations in terms of theta functions 



a 


= 


[U]{U, X}{U,Y}{U,Z} 


bi 




= [X]{X,Y}{X,Z}{U,X} 


cii 


= 


[U,X]{U,Y}{U,Z} 






= [X,Y}{X,Z}{U,X} 


a\ 


= 


[U,Y]{U,Z}{U,X} 


b> 2 




= [X,Z}{X,Y}{U,X} 


a'l' 


= af = 


[U,Z]{U,X}{U,Y} 


bo 




= [X,U]{X,Y}{X,Z} 


a 2 


= 


[U,X + Y}{U,Z} 


bs 




= [X,Y + Z]{U,X} 


4 


T 

= «2 = 


[U,Y + Z]{U,X} 


b'(' = 


b\ 


= [X,U + Z]{X,Y} 


A 


= 


[U,X + Z}{U,Y} 


b\ 




= [X,U + Y]{X,Z} 


a 3 


= 


[U,X + Y + Z] 


bt= 


b T 

u 2 


= [X,U + Y + Z] 


c 2 




[Y]{Y,Z}{U,Y} {X,Y} 






= [Z]{U,Z}{X,Z}{Y,Z} 


c\ 




[Y,U]{Y,Z}{X,Y} 


4" = 


df 


= [Z,U]{X,Z}{Y,Z} 


c 3 




[Y,Z]{U,Y} {X,Y} 


d' 2 




= [Z,X]{Y,Z}{U,Z} 


Cl 




[Y,X] {Y,Z}{U,Y} 


d 2 




= [Z,Y]{U,Z}{X,Z} 


Co 




[Y,U + X]{Y,Z} 


d'( = 


d\ 


= [Z,U + X]{Y,Z} 


Jtt 


= C% = 


[Y,Z + U]{X,Y} 


d\ 




= [Z,X + Y]{U,Z} 






[Y,X + Z]{U,Y} 


d\ 




= [Z,U + Y]{X,Z} 




= c\ = 


[Y,Z + U + X] 


d 




= [Z,U + X + Y] 



Table 1. The 32 line section parameters appearing in eas. 1)19(1 . expressed in terms of 
theta functions and prime factor ratios, using the abbreviations l)26[l. Observe that in 
the prime factor brackets { , } the points come always in the order U, X, Y. Z (without 
and with primes). 



For writing the MTE in our new parameterization and to check (jl8)l and (fTT^) . we 
consider three more spaces V = C N , corresponding to the indices 4, 5, 6 . In Figure 121 
the first three spaces were located at the intersection points of the lines X, Y, Z . To 
include the other three spaces, consider the "quadrangle" formed by four lines shown 
in Figure El Corresponding to the new line U we introduce another pair of points 

u',u er g . 
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U 



Figure 3. Quadrangle in the auxiliary plane formed by the directed intersection 
lines of four oriented lattice planes. The six spaces V in which the MTE operates are 
considered to be located at the six intersection points. 



Looking at Figure El we see that, instead of labeling the spaces by the vertices 
1, . . . , 6 of the quadrangle, we can as well label them by the pair of lines which inter- 
sect in these vertices, so identifying 

1 ~ YZ; 2 ~ XZ; 3 ~ XY ; 4 ~ UZ ; 5 ~ UY ; 6 ~ UX . (36) 

Note that the ordering of the lines is important for the identification (J3T?|) : we shall chose 
the anti-clockwise orientation in Figure HU not the mirror reflected clockwise one. 

Next we assume that we can write the "coupling constants" Kj all in the form (|33|) . 
Then from (|3l)|) it is suggestive to build e.g. k± from the points Y', Y, Z', Z only, etc. 
and put (factors q 1 ^ 2 are inserted to produce correct signs when forming N-th powers 
for jTH)): 

1/2 v'>^} 1/2 {^''^} 1/2"P^')^} 

H {U,Z} ' H {U,Y} ' H {U,X} y 1 

Now we consider (|3*T|) and (|3*7jl and assume that the line-section parameters ao, cti, ■ ■ ■ 
and do, d%, . . . follow the same scheme as postulated for bo, ... , cq, . . . above after 
a, ~ [C/,...], di ~ [Z, . . .]. So equations (|TS|) lead us to express all line-section 
parameters in terms of theta-functions as shown in Table ^ We make ample use of the 
short-hand notations (J26j) . 

Apart from Q which always comes with v, we use eight arbitrary points X', X, Y', Y, 
Z', Z, U', U (zT g . The Kj may as well be written in terms of the bared brackets using 
{A, B'}{A, B} = {A,B}{A',B} . From Table □ we see that the a { don't depend on 
U', the bi not on X' etc. 

Now, using the results of Tableland (J37j) . we shall re- write all the Hirota equations 
(J19|) in terms of theta functions on T g and prime form cross ratios. Not very surprisingly 
in view of [UJ [T71 [THj, it turns out that these all have the form of the double- Fay 
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identity. Also, as expected from Figure E3 and the meaning of the mappings as moving 
lines within the quadrangle, the 24 equations follow from each other by a sequence of 
simple substitutions. Just inserting from Tableland (JBZJ), the first three equations of 

(US)) become, using the notation (|2*7j) (recall that these are the equations (fT2"j) defining 

( f) 

the functional mapping TZ^): 



F(v; X,Y',Y,Z',Z) 



F(v + I$';Y',X',X,Z f ,Z) 



({U, X}{U, Y}{U, Z}{X, Y}{X, Z}) 
E(X, Y) E{X, Z) E{Y, Z) 

{{UjC}{U^}{U^}) 
E(X, Z) E(X, Y') E(Y',Z) 



N 



N 



F(v; Z,Y',Y,X',X) 



({U, X}{U, Y}{U, Z}{X, Z}{Y, Z}) 



N 



0. (3* 



E(X, Y) E(X, Z) E(Y, Z) 

The dependence on U', U appears only in the factors on the right, not in the F . 
Assuming generic points U', U, . . . we conclude that the F must vanish and we combine 
the essential terms of into 

/F(v;X,Y',Y,Z',Z) \ 



3(v; X',X,Y',Y,Z',Z) 



def 



(39) 



/ 



F(v + I^'; Y',X',X,Z',Z) 
\F(v; Z,Y',Y,X',X) 
Then the 24 Hirota equations (jl9j) which describe the functional mappings take the form 
3 (v; X', X, Y', Y, Z', Z) = ; 3 (v + l u J ■ X', X, Y>, Y, Z', Z) = ; 

ff(v + lf; U'.C/.n^'^) = 0; ff(v; C/', Z7, F', Y, Z', Z) = 0; 



0; 



(40) 



^ (v; C/', U, X', X, Z', Z) = ; (v + 1^'; U', U, X', X, Z', Z) 

y(v + i|' ; f/',f/,x',x,r,y) = o ; y(v ; t/'.^x',^,^,^) = o. 

The equations in the left column of (j4*Uj) precisely correspond to the first twelve equations 
of (fTTIjl . The last three equations of (fTT)|) are combined into the top equation of the right 
column of (|4~0|). As already mentioned with (|2~Uj) . equations pUj) together contain the 
functional TE. 



2.4- Theta-parameterization of the simple modified tetrahedron equation 

Fianally, we use the parameterization (}3*Tj) to re-write the MTE. From (j!7|) with (fTHj) 
we find that the functional mapping just produces a permutation of the four pairs of 
points X, X', . . . ,U,U f G T g , together with shifts in the vector v. Of course, the result 
corresponds to (JJ0J). Explicitly it is: 

Theorem 1 The simple modified tetrahedron equation may be parameterized in the 
terms ofY g , v e Jac(r 9 ) and four pairs X' ', X , Y',Y, Z',Z, U',U G T g by the definition 
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W>W>W as follows: 

R (1) = R(v; X', X; Y', Y; Z', Z) ; R (5) = R(v + l£f; X', X; Y', Y; Z', Z) ; 

R (2) = R(v + if; U', U; Y', Y; Z', Z) ; R (6) = R(v; U', U; Y', Y; Z', Z) ; 

(41) 

R (3) = R(v; U', U] X', X; Z', Z) ; R (7) = R(v + Iy ; U', U; X' , X; Z', Z) ; 

R (4) = R(v + if; U', U; X', X; Y', Y) ; R (8) = R(v; U', U; X', X; Y', Y) . 

Proof: Each R^) is determined by its three Fermat points n . From 

[Oj these points are known in terms of the line-section parameters, see ()18j) . Inserting 
the theta-function expressions for the line-sections from Tabled into (fT%j) one finds that 
the xf' for j = 2, . . . , 8 are obtained from those for j = 1 , equations (jSIJ), by the 
substitutions seen in (jJTJ). □ 

Using the correspondence between the labels 1, . . . , 6 and the line labels U, X, Y, Z, 
given m flUJ), R« = R123 may also be labeled as R XYZ etc., and we write the MTE 

as 

R XYZ (v) R UYZ (v + If) R uxz (v) R UXY (v + If) 

= p R UXY (v) R uxz (v + l Y y) R UYZ (v) R XYZ (v + $g). (42) 

This notation also indicates directly the three pairs of points on the algebraic curve 
which parameterize the matrices R^ in f)41j) . 

In jU] we had shown that using simple re-scalings, out of the 24 line-section 
parameters listed in Table E] and the 6 parameters K\, . . . , , only eight parameters 
are independent. Here we have eight points on T g which can be chosen freely. In 
addition, 16 phases from taking the N-th roots can be chosen freely. In terms of the 
line-section parameters, the choice of the independent phases is the same as the one 
explained in 0. 



3. The fused vertex weight 9^ 

3.1. Open Ni x N 2 x N 3 box 

The natural graphical interpretation of the R-matrix is a three dimensional vertex, i.e. 
the intersection of three planes in 3-D space. The six indices aj, a'j are associated to 
the edges of the vertex, recall Figure [T] 

The next step is the consideration of the intersection of three sets of N±, N 2 and 
N 3 parallel planes. This produces a finite open cubic lattice of the size Ni x N 2 x N 3 . 
We call the corresponding vertex object 9^. It is the result of the fusion of elementary 
R-matrices. 

The lattice is defined as in ((TJ), but now rij — 0, . . . Nj — 1 . Let 9ti 2 3 be the matrix 
associated with the open cube (more precisely, the open parallelepiped): 

9^123 = (<?1, <?2, <?3| 91 \a[, a' 2 , <? 3 ) = / ] \ \ (g"l,n; 0"2,n+e 2 ) °"3,n|R-nkl,n+ei ; &2,m Q"3,n+e 3 ) • 

M " (43) 
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Here the six external multi-spin variables (i.e. the indices of the matrix 9^123) are 
associated with the six surfaces of the cube: 

C?l = {°"l:0,n 2 ,n 3 }, <?2 = {^im ,7V 2 ,n 3 } , <?3 = {^3:n u n 2 ,o} , Uj — 0, Nj — 1 (44) 

and 

V\ = {Cl-.m^ns}, &2 — {o^mAns}' ^3 = {^m.na.ATg}, % = 0, Nj — 1 , (45) 

and the summation in (j4Hj) is taken with respect to all internal indices o~j >n . Anticipating 
what will be needed in (J62J) in order to prove that the fused weights satisfy a MTE of 
the same form as we had in (|4*2|) . we use a reversed numbering for o~i and a' 2 , so that 
the "initial" external indices are <j\ = 0, o~ 2 = N 2 , 03 = 0. This reversed numbering 
in the second space is dictated also by our choice of the line orientations in the lattice 
and, as consequence of this, in the auxiliary plane (see Figure HJ). 




Figure 4. Top: 3-dimensional view of the oriented cube N\ = N 2 — N 3 = 2 (heavy 
lines) which is formed by six planes (indicated by dashed lines) . Bottom: the horizontal 
auxiliary plane with the three pairs of lines arising from the intersection of the three 
pairs of planes with the auxiliary plane. This is a generalization of Figures and [21 If 
we consider e.g. the point (0, 1,0) to correspond to the point A of Figure QJ then the 
inner triangle in the auxiliary plane corresponds to the left shaded triangle of Figure ^ 
and to the left part of Figure [21 So the initial external spin (Weyl) variables (|44|l can 
be considered to sit at the three times four intersection points of the auxiliary plane. 
In order to get a similar picture for the final external variables we have to place the 
auxiliary plane above the cube. 
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In our next step we want to parameterize all R n in (|43|) such that the fused 
weight D\ again satisfies a Modified Tetrahedron Equation. We shall show that using a 
theta-function parameterization this is possible and the ttijk obtained will depend on 
Q(Ni + N 2 + N 3 ) free parameters. 

We use again the generic algebraic curve T g , and one vector v G C 9 . As in (J3*U|) 
each R n will depend on three pairs of points on T g , and to each R n we assign different 
three pairs: 

X' ni , X ni , Y' n2 , Y n2 ,Z' n3 ,Z n3 , rij — 0, Nj — 1. (46) 

However, the argument v will be shifted for each R n by an amount I n which depends 
on the points assigned to "previous" neighbors: We define 

r£*> = R ( v + l n; x^,X ni -X 2 ,Y n2 ;Z' n3 ,Zj , n, (I V, - 1. (47) 

where 

m— 1 ri2— 1 "3 — 1 

in = E X S; + E + E x £: • ( 48 ) 

mi=0 m2=0 rri3=0 

Now (J4*3"j) and (|4"7)l define the matrix function 

9Wv) = ^(v;X',X;y',F;Z',Z) (49) 
where X', X, Y', Y, Z', Z stand for the ordered lists of divisors, 

X = (X ,X 1 ,...,X^ i _ 1 ), X' = (X^,X;,...,X^ i _ 1 ), Y = {Y ,Y U ...), etc. (50) 

As to the index structure, recall ©, 

m 123 e En.d(V NaN * ®V NlN * ®V NlNa ) , (51) 

where in the same way as before we enumerate the number of V 3 k in the tensor 
product (e.g (|4*H|l are the matrix elements of ^123). In Figure0we show the intersection 
lines of the planes of a Ni x N 2 x N 3 cube which appear in an auxiliary plane (as in 
the bottom part of Figure Ej), which intersects the "initial" edges corresponding to (jUJ). 
On the section the jVj + N 2 + N 3 planes become lines, and the edges of the cubic lattice 
become the intersection points of N\N 2 + N 2 N 3 + N\N 3 lines in the auxiliary plane. 
The intersection points are gathered into three sets Vi = V^ 2 ^ 3 etc., and the index of 
Vj is the number of the corresponding ! in the tensor product in ()51|) . Figure El is 
helpful to arrange the numbering in (|4*Hj) and the correct assignment of X' , X ni , etc. 

3.2. The Modified Tetrahedron equation for the fused weights 

For writing the MTE, apart from the three pairs of sets X',X; Y',Y; Z',Z of (j4*9*|) . 
fl5()(l we need a fourth pair 

u = (u , u u ...,t/iVo-i); u' = {u' , u[, ...,^ _ 1 ). 

Applying the definition (|43|h in addition to ()49|) we construct the matrices 

^145 (v) = 5H(v; [/', [/; Y', Y; Z', Z) ; ^ 246 (v) = *H(v; C/', ?7; X', X; Z', Z); 




Figure 5. Ordering of the indices of the matrix EH123, shown for the case N% = N2 = 
N3 = 3 by drawing the auxiliary triangle in the auxiliary plane, as in the bottom part 
of Figure 



Wv) =9t(v;tf' \U;X',X;Y',Y). 
Their index structure is denned by 

Vl = V N 2 N :i ) y 2 = yJVstfi ? ^ = y^iVa ? 

v 4 = V^^ 3 , V 5 = V N ° N2 , V 6 = V NoNl . ^ } 

so that e.g. 9*145 is acting in a space of dimension 7V^ 2 3+ ^°^ 3+Ar °^ 2 . For 9*145 in 
analogy to the definitions P7 j) .()48| l one uses 

= R(v + I n ; U^U no ; Y n2 ; Z' ns , Z ns ) (53) 

with 

no— 1 n 2 — 1 "3 — 1 

: » = E lU ut° + E 1^ + E l tl > ( 54 ) 

mo=0 m,2=0 ni3=0 

similarly for 9* 2 46 an d 9*3 56 . 

Theorem 2 T/ie matrices 9* defined in UJfy obey the modified tetrahedron equation 

9*123 (v) 9*145 (v + Ix) 9^ 2 46 (V) 9* 356 (v + Iz) = 

(55) 

p 9*356 (V) 9* 2 46 (V + Iy) 9*i 45 (v) 9*i 23 (v + Ip) 

where 

T _ V^JVo-l T^O . T — ST Nl ~ l T X ™1 ■ 

L U — 2^n =0 L U no ) A ^ — Z^n 1= A X rei > 

(56) 

T V^A^2-1 T^"2 . T V~^JV3— 1 T^"3 
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Proof: The main content of this theorem is the appearance of the specific set of 
shifts f!48|) . (jHSI)- For the proof it is convenient to introduce some compact notations. 
Instead of using the number labels for the 9* we shall use the labels U, X, Y, Z just as 
these were introduced in (jHBj) and (pI2|) for the single vertex matrices R. So, for the box 
9t-matrix and its sets of divisors we write, 

9*123 (v) => 9* XYZ (v); <K 145 (v) => <K UYZ (v); etc. (57) 
In this short notation , formulas (pIHj) . (jTfj) imply the definition 

fH xYz ( V ) = n n 1 1 R^^^Kv+in) (58) 

711=01^1-1 ^2=^2-110 713=01^3-1 

where we use ordered products 

JJ fni = fofl • • • fiVi-l , Yl f»a = fjv 2 -i • • • fifo ■ (59) 
ni=0TJVi-l 712=^2-1^0 

For the triple (X, Y, Z) , I n is given by Each R Xn i Yn 2 z ™3 acts non-trivially in only 
three of all the spaces (|5l?J) and the ordering is relevant just for neighboring indices X ni 
or Y n2 or Z n . A . The analogous modifications required to get the other matrices 9*145, etc. 
should be evident. 

Now we turn to the proof of the theorem which will be by mathematical induction. 
The theorem claims the validity of the MTE for arbitrary N , N\, N 2 , N 3 . For the initial 
point Nq = Ni = N2 = N3 = 1 the MTE holds because it is just Then to prove 

the theorem, we reduce the MTE f!55|) for some Nj to MTEs with iVj < Nj. Thus one 
has four similar steps of the induction. Here we consider the induction step for the 
X-direction, the other steps follow analogously. 

We split the list X into two sublists of length and N[ 2) = JVi - Njp: 

X« = (X ,X 1 ,...,X N m_ 1 ) , X^ = (X N m,...,X Nl ^) , (60) 

jV (1) -l 

so that l x w = z2ni=o ^-x n an d I X (2) = Ix — I^C 1 ) ■ According to this splitting and due 
to the definition (J58j) 

v) ^ x(2)YZ (v + I x( i)), 
uxWz 



9* XYZ (v) 


= 9t x(1)YZ ( 


<K uxz (v) 


= 9* ux(2)z ( 


<K UXY (v) 


_ ^UX( 2 >Y 



v + I x( i))9* uxwz (v), (61) 

[v + l xm ) 9* ux(1)Y (v). 

The meaning of notations X^ 1 ) and X^ 2 ^ used in (161)1 should be evident from ([58)1 . Observe 
that because of the reverse numbering with respect to the middle superscript of 9t in 
()58jl . the factors in the latter two equations appear in reverse order. Since 9t UYZ (v) 
contains no X, neither as subscript nor in the argument, it will not be split. However, 
we have to put 

9* UYZ (v + I x ) = 9* UYZ (v + I x(1) +I x(2) ). 
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Now substituting (JoTj) to the LHS of (fo3|) written in our new superscript notation, and 
abbreviating V! = v + I X (V) , we get 

<K XYZ (v) <K UYZ (v + l x ) <K uxz (v) <K UXY (v + l z ) 

(v + l xw ) ^ UYZ (v + I xm + I x(2) ) lH ux(2)z (v + I xm ) X 



^XCX)YZ (v) ^>YZ, 



x 9* ux(1)z (v) 9^ ux(2)Y (v + l xw + l z ) !H ux(1)Y (v + l z ) 



^ x(1)YZ (v) 



X( 2 >YZ, 



Vl ) <K UYZ ( Vl + I x(2) ) 9^' z ( Vl ) Y ( Vl + I z 



UX( 2 )Z/ 



>UX( 2 )Y, 



^ ux(2)Y ( Vl ) m 



x$n ux(1)z (v)^ ux(1)Y (v + i z ; 



ux(2)z (v 1+ I y ) 9^ UYZ ( Vl ) ^ x(2)YZ ( Vl + i c 



x 



X 



x$H uxll)z (v)lH ux(1,Y (v + I z ; 



^ ux(2)Y ( Vl )^ ux(2)z (v 1 + I y )x 



X 



^x(x)yz (v) ^uyz (v + ^uxwz^ ^uxW Y(v + ^ 



^ ux(2)Y ( Vl ) ^ ux(2)z (v 1+ I y )x 



X 



^uxWy (v) w uxwz (v + Iy) ^uyz (v) j^xWyz^ + T , 

^UX(2)Y (v + g^UXWY^ ^UX(2)Z (v + ^ + j y) ^UX^Z 



^ x(2)YZ (v 1+ Ia) 



x^ UYZ rv)9l x(1)YZ (v + Ir/)5H x(2)YZ i 



V + Iy) X 

V + I x( i) + Lj) 



UXY, 



(v) <K uxz (v + I y ) <K UYZ (v) 9T YZ (v + l v ) . (62) 

From the third to the fourth line of ()62j) within the inserted brackets we used the MTE 
for the smaller set (U, X@\ Y, Z). In order to be able to isolate the terms containing 
the order of factors in the second line of (jfilj) . which was used in the first step, 
is crucial. Going from the fifth to the sixth line in the brackets we used the MTE for 
(U, Y, Z) . In the other steps of (}62|) we just commuted or combined various terms. 

The last step is made possible by the "reverse" order of factors in the last line of ([61)1 . 

□ 



4. Special cases: Solving the tetrahedron equation 

4-1. Compact algebraic curve 

Let now N = N± = N 2 = N 3 = M . Starting from the generic parameterization of the 
MTE the usual tetrahedron equation is obtained if 

9t(v) = 9t(v + 1) and p = 1 . (63) 

This is the case when 

l v = I x = Iy = l z = on Jac(r 9 ) , (64) 
and the ratios of ^-functions (125)1 are periodical. 
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According to Abel's theorem, (J64j) means that there are four meromorphic functions 
u, x, y, z on T g with the divisors 

M-1 M-1 M-1 M-1 

( U )=Y, U no- U n , (X)=^A; 1 -A ni , (y)=J2Y; 2 -Y n2 , (z) = £ K z ~Z«*. 

710=0 711=0 712=0 713=0 

(65) 

As it is well known (see e.g. Theorem 10-23 of |2()j). conditions ()64)) are a strong 
restriction for the type of T g : T g is the algebraic curve given by the polynomial equation 

M 

P(x,y) d = f x a y b Pa>b = 0. (66) 

a,b=0 

The choice of any pair of x, y, z, u produces an equivalent polynomial equation. The 
form of the polynomial P(x, y) provides the restriction for the genus, 

g < (M-1) 2 . (67) 

Thus we come to 

Theorem 3 Let T g be the compact algebraic curve defined by polynomial equation A 6' 6)) . 
Let four sets of U' no ,U no , X' ni ,X ni , Y' n2 , Y n2 and Z' m ,Z m , n k = 0,...,M- I, are the 
divisors of four meromorphic functions u,x,y,z on T g . Then the tetrahedron equation 
is satisfied 

9ii23(x, y 7 z) 9ti45(u,V,z) ^24g(u,x,z) <R 3m (u,x,y) = 

(68) 

91 3 5 6 (u, x, y) 9i24 6 (u, x, z) 9*i 45 (u,y,z) 9*123 (x, y,z) , 
where four matrices are the same matrix function of different arguments, 

9K{x,y,z) = 9*(v; X\ X; Y', Y; Z', Z) etc. (69) 
defined via HJW, HJ^), HTW an d ( E^ )- 

According to the conventional terminology, one may say that u, x, y, z are the spectral 
parameters, the moduli of T g are the moduli of the tetrahedron equation, and vector v 
is a kind of deformation parameter. 

Theorem El may be also be proved differently, without mentioning the simple MTE. 
In this alternative approach one considers the auxiliary linear problem for the whole box 
and the corresponding mappings. See e.g. [2H 122] for the description of this method 
and [23J for the parameterization of the classical equations of motion. Remarkably, in 
this alternative way the spectral curve (f6T)|) appears naturally from the linear problem. 

4-2. Simple tetrahedron equation for ZBB-case recovered 

In section 1 we discussed the ZBB- model and its R-matrix (JJJ). We now show how our 
scheme contains this case. Formally ZBB's tetrahedron equation corresponds to 
with M = 1. It gives g = 0, i.e. the spectral curve is the sphere with 

E(X ' Y) = 73TW' e() ^' (70) 
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Here the formal theta-function has no argument since the Jacobian is O-dimensional. 

Conditions ()64j) therefore are out of use, and the parameterization ([30 )1 .([3i p contains 

(X — Z)(X' — Z') 

the N-th roots of the cross-ratios like -, -. One may show that the 

(X' - Z){X - Z') 

number of independent cross-ratios is the number of variables X, X', ... minus three. 
Therefore the single R-matrix contains 6 — 3 = 3 independent complex parameters (as 
it should be), and the simple tetrahedron equation contains 8 — 3 = 5 independent 
complex parameters (again as it should be). It means that the tetrahedral condition, 
which appears in Zamolodchikov's parameterization of R in the terms of spherical 
triangles, is taken into account automatically. Moreover, the parameterization with 
the help of the cross-ratios takes automatically into account the geometric structure of 
any set of the planes in three dimensional Euclidean space. The parameterization of 
the inhomogeneous Zamolodchikov-Bazhanov-Baxter model in the terms of cross-ratios 
corresponding to the g = limit of (p")J) . (j3Tf . (pTTj) has already been used in [2I|. 



4-3. Chessboard model 

Previously derived "chessboard models" of the lattice statistical mechanics based on 
the modified tetrahedron equation ^UJ are of course related to our considerations. 
The term "chessboard" appeared as the visual interpretation of the cubic lattice with 
M = 2 to be homogeneous. It means that the cubic lattice consists of eight different 
types of vertices (i.e. eight different types of the Boltzmann weights) - a kind of three 
dimensional analogue of the chess board with eight different colors of the cells. 

The models described in [TUl HI] are at the first the so-called IRC-type models, 
but with the help of vertex-IRC correspondence [Hj one may construct their vertex 
reformulation. Thus the model implicitly described in ^T] is equivalent to M = 2, 
g — 1 of our scheme. For g = 1 the curve and its Jacobian are isomorphic, so that 
without loss of generality one may chose 

I* =X-Y G C/Z + Zr. (71) 

Further, one may use 9\ 

oo 

e(v) = O^v.r) = e CTr( " +1/2)2+2 ^ +1/2)( " +1/ ^ (72) 

n=— oo 

as basic theta-function, and E(X,Y) ~ 6i(X — Y) as the prime form. This formulas 
simplify the definitions ()25)1 . Periodicity conditions (|65|) may be chosen 

Xq — Xo + X ; x — Xi = Yq — Y + Y' : — Yi 

= Z' -Z + Z; - Z x = % - U + U' x - Ui = 1 . (73) 

Note, the 1 in the right hand side of (JTSJ) is equivalent to r (due to the Jacobi transform), 
while instead of 1 in the right hand side of ()73|) gives a trivial model. 

The model explicitly described before in |l(Jj corresponds to M — 2, g — 1, X — X^, 
X^ = x;, Y = Y x , Yq = Y[, Z = Z 1; Z = Z[ etc. with the condition (J73J). This choice 
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leads to the identification of the parameters in 

-R-n -R-n+ei+e2 ^n+ei+ea -R-n+e2+e3 > C^) 

so the cells of this three dimensional chess board have only two "colors" . 

Note that the vertex-IRC duality is not exact because it changes the boundary 
conditions. 



5. Conclusions 



We considered a large class of integrable 3-D lattice models which have Weyl variables at 
N-th root of unity as dynamic variables. We have shown how the Boltzmann weights can 
be conveniently parameterized in terms of N-th roots of theta functions on a Jacobian 
of a compact Riemann surface. The Fermat relations of the points determining the 
Boltzmann weights are simple Fay identities and the classical equations determining 
the scalar parameters are a consequence of a double Fay identity. In the modified 
tetrahedron equation we have four pairs of arbitrary points on the Riemann surface in 
simple permuted combinations. 

This parameterization allows a compact formulation of the rules to form fused 
Boltzmann weig hts m e End C 3NM which are the partition functions of open boxes 
of arbitrary size. The d\ obey the modified tetrahedron equation and are again 
parameterized terms on iV-th roots of theta-functions on the Jacobian of a genus 
g — (M — l) 2 compact Riemann surface T g . The spectral parameters of the vertex 
weight !tH are three meromorphic functions on the spectral curve T g . For the case 
that the Jacobi transforms become trivial the 9i obey the simple tetrahedron equation. 
The Zamolodchikov-Baxter-Bazhanov model and the Chessboard model are obtained as 
special cases. 

So, the scheme discussed here contains and generalizes many known 3-D integrable 
models. The hope is that the framework is now sufficiently general to contain physically 
interesting models with a non-trivial phase structure. However, to get information on 
partition functions, either analytically or approximately, is still a very difficult open 
problem. There is no way known to generalize Baxter's special method by which 
he obtained the partition function of the ZBB-model. 
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